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Abstract 


Fuel cell applications become increasingly attractive. Therefore, comprehensive models, simulation software and analysis tools are required to 
characterize fuel cell behavior. This paper introduces a model capable of representing the static and dynamic (dynastatic) behavior of PEM fuel cells. 
The model approximates the theoretical current-voltage description with an equivalent electric circuit. The parameter values of the model are found 
by analyzing small-signal measurement results from electrochemical impedance spectroscopy, together with the steady-state electric characteristic 
and large-signal behavior from step-response measurements. Because the modeling approach is completely based on measurements on the fuel cell 
electric output terminals, electrochemical aspects like temperature dependency, charge double layer and adsorption effects are implicitly covered. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Fuel cell applications become increasingly attractive, because 
they are efficient, clean, and, in the case of a polymer electrolyte 
membrane (PEM) fuel cell, light weighted electricity genera- 
tors [1-9]. Because the interest in fuel cells grows, the urge for 
accurate simulation models grows as well. 

The focus of the analysis and modeling of fuel cells has 
always been on steady-state behavior. However, fuel cells are 
low voltage, high current electricity generators, which is not 
practical in common applications. By using power electronics 
to convert the low voltage to a more efficient high level, one will 
introduce step currents, and other dynamics in the output of the 
fuel cell, which can certainly not be modeled by one of the two 
known steady-state models presented in the literature [2-9]. 

In the literature a transient model based on operation temper- 
ature changes was already described [4]. A major disadvantage 
of that model is that it requires temperature measurements. 
Commercial available fuel cells stacks are not always equipped 
with (enough) temperature sensors. However, it is possible to 
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characterize the fuel cell dynamics also by applying step load 
changes [9,10] or small current variations at the output terminals 
[9-12]. Unfortunately, no qualitative nor quantitative compari- 
son between both approaches can be found in the literature. 

This paper presents three measurement methods (Section 
2) and accompanying test results (Section 3), namely electro- 
chemical impedance spectroscopy for small-signal behavior, 
steady-state measurements to determine the current-voltage 
characteristic and step-response to characterize large-signal 
behavior. A unique equivalent circuit model based on these three 
measurement sets is introduced in Section 3, and is validated 
in Section 4. This model is capable of simulating static and 
dynamic fuel cell behavior (a so-called “dynastatic” model) and, 
because this model is completely based on electrical measure- 
ments of the output terminals, it implicitly covers temperature 
effects, charge double layer and adsorption dynamics. 


2. Method and materials 


The fuel cell voltage, Vac, can be expressed as a function of 
the fuel cell current Jpc [8-10] 


I . I 
Vec = Eoy — A In (=) R™ Ipc + B In (1 = že) (1) 
0 1 
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Fig. 1. Schematically drawn PEM fuel cell, with (a) the anode separator plate, 
(b) the anode with a carbon-based diffusion layer and platinum catalyst, (c) is 
the polymer electrolyte, (d) is the cathode with a carbon-based diffusion layer 
and platinum catalyst and (e) the cathode separator plate. 


in which Eocv is the open circuit voltage of the electrochemi- 
cal reactions. Furthermore, three loss terms are included in (1). 
The first is the activation loss, which models the slowness of the 
reaction at the surface of the electrodes and is characterized by 
parameters A, the activation coefficient, and Jo, the exchange cur- 
rent. The second loss term is the internal resistance R'™, which 
models the ohmic losses in the conductors and the membrane. 
The third term is the concentration loss, which models the limit 
of the fuel flow to the active area of the membrane. The concen- 
tration loss is characterized with B, the concentration coefficient, 
and J, the limiting current. 

From the mechanical construction of a PEM fuel cell, shown 
in Fig. 1, one can imagine that two ideally polarized metal 
plates separated by a thin plastic sheet act as a capacitor. In 
case of the PEM fuel cell, the large-signal dynamics of the sys- 
tem can be modeled using a capacitor in combination with a 
parallel resistor and a series resistor [8]. However, a fuel cell 
is an electrochemical device, and the referred capacitance is 
caused by a charge buildup between two different materials such 
as the electrode and electrolyte. As is shown in Fig. 1, a fuel 
cell contains two electrode—electrolyte interfaces, namely the 
anode—membrane and the cathode—membrane interface. There- 
fore, it may be expected that a fuel cell suffers from two charge 
double layer capacitances, resulting in two time constants [13]. 

Moreover, the platinum catalyst used to speed up the anode 
and cathode reactions, may result in an adsorption time constant 
[13-15,10]. This can be explained by the fact that hydrogen 
adsorbs to the platinum catalyst surface, as a result of hydrogen 
ion discharge. Hydrogen has the property to transfer directly 
almost completely its electron charge to the catalyst [16], result- 
ing in an immediate voltage response, while the current of the 
system, e.g. the transport of protons and electrons, is not directly 
influenced, because the proton transport through the membrane 
is slow. This phenomenon creates a phase shift between current 
and voltage, hence, an additional time constant [9,10]. 

We performed three types of measurements on the fuel cell. 
First, impedance spectroscopy finds the small-signal response 
of the fuel cell in certain DC set-points. Second, steady-state 
measurements are done to determine the voltage—current char- 
acteristic and the maximum power rating of the fuel cell. Last, 
step-responses reveal the large-signal behavior of the fuel cell. 


The fuel cell used for the measurements is 1 of the 12 car- 
tridges from a commercially available 500 W fuel cell system 
[17,18]. The cartridge contains four membranes stacked in pairs. 
For the experiments, one cartridge is used and the two pairs of 
membranes are connected in series to increase the output volt- 
age and limit the output current. The cartridge is mounted on a 
aluminium frame with a fan, which provides the oxygen to the 
cathode of the fuel cell. The fuel cell in this test-setup delivers a 
maximum output power of only 7 W due to insufficient airflow 
supply. 

To perform the measurements, the fuel cell is coupled to 
a linear regulator, which is controlled by a digital signal pro- 
cessor [10,9]. The voltages over the drain and source of the 
regulator, and the current of the fuel cell are recorded by 12 bit 
non-multiplexed AD converters. The accuracy of the digitally 
measured voltages and currents is limited by the bit-error of the 
used converters. 


3. Results 
3.1. Electrochemical impedance spectroscopy 


Electrochemical impedance spectroscopy (EIS) is a method 
to characterize electrical properties of materials and their 
interfaces with electronically conducting electrodes [13]. The 
measurement procedure is to excite the output terminals of a fuel 
cell at a fixed DC set-point with superimposed small amplitude 
sinusoidal signal. The complex impedance can be interpreted as 
an electrical circuit containing resistors, capacitors and induc- 
tors. Each time-constant is a result of a chemical, physical or 
electrochemical phenomenon. 

The impedance spectroscopy measurement results carried out 
on the fuel cell stack are shown in Fig. 2. The frequency ranges 
from 0.01 Hz (right side of the semi-circles) to 45 Hz (left side 
of the semi-circles), and the DC-currents were set to 1, 2, 3, 4, 
5and 6A. 

In Fig. 2, the averaged data-points obtained from the measure- 
ments are shown. The test is repeated 10 times for all frequencies 
and for all DC points. The data is converted to the frequency- 
domain using a fast Fourier transform. 


Z" (Q) 


Fig. 2. Impedance plots of the fuel cell in a range from 1 to 6 A. Z’ and Z” are 
the measured real and imaginary impedance components, respectively. 
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Fig. 3. Third-order ladder network for fitting impedance spectroscopy data. 


The measured impedance show a large capacitive semi-circle 
in the frequency range of 0.1-45 Hz, which is caused by charge 
double layer effects. Furthermore, inductive behavior is seen in 
the low-frequency range, which is the right side of the impedance 
semi-circles. References in the literature indicate that this behav- 
ior is probably caused by the adsorption of hydrogen to the 
platinum surface [14,15,10]. 

The measurement results of the electrochemical impedance 
spectroscopy show inductive and capacitive behavior, which 
implicates that the equivalent electrical impedance should con- 
tain at least one capacitor and one inductor. The proposed 
equivalent impedance circuit, with two double layer capacitors 
CpL and CR and one adsorption inductor La, is shown in Fig. 3. 
A previous publication indicates that the fuel cell stack used is 
best represented with a third-order equivalent circuit model [9]. 

The the non-linear least-squares fit provides an accurate 
description of both the amplitude and the phase-shift. The values 
of the equivalent model are presented in Table 1. 

The resistors found with EIS are graphically represented in 
Fig. 4. It is clear that the bulk resistor Roo has a constant char- 
acter, and can therefore be approximated as a constant with 
value Rog = 0.272 Q. The adsorption resistor Ra tends to slowly 
decrease with increasing current set points. Therefore, this resis- 
tor is approximated using 


RaUpc) = RÌ — æa Irc, (2) 


where RÀ, = 23.5m Q represents the initial value and aq = 
2.62m Q AT! is the linear coefficient. 

The two double layer resistors (Rp and Rp) are showing 
non-linear behavior dependent on the fuel cell current. From 
Joule’s law, it follows that the temperature dependency of a con- 
ductor is proportional to the square of the current through the 
conductor, that is, 


T œ žc- (3) 


Table 1 


Therefore, in first approximation a function representing the 
resistors is proposed as being proportional to the square of the 
fuel cell output current. The implementation of the resistor func- 
tions is 


RoiUgc) = R + RY, eebe, (4) 
and 
Rr(Ipc) = RX + Roe Rc, (5) 


The settling value for the first double layer resistor is R}, = 
75.2m Q, which represents the resistance for large output cur- 
rents (Ipc). Furthermore, the initial value Rey = 172m Q and 
the exponential rate coefficient ap}, = 0.342 A~?. The second 
double layer resistance is characterized with RR = 18.4m Q, 
RÌ, = 97.0m Q and a rate coefficient a = 0.244 A~*. The fit 
is presented in Fig. 4(a). 

However, from the measurement data presented in Fig. 4, the 
exact characteristic of the double layer resistors cannot be deter- 
mined in the low current region. This is due to the limitations 
of the measurement-setup to measure the injected small current 
variations so that the obtained voltage variations are undistort. 

A second approach assumes that the temperature changes are 
small. In that case, the current-temperature relationship (3) can 
be linearized, resulting in 


T x lfc. (6) 


Using this linear relationship, the resistor functions can be 
represented with 


RpL(rc) = RỌ + RO e LHe, (7) 
RrUgc) = RY + Rae RC, (8) 


The values for the first double layer resistor are Rhy = 67.8 m 
Q, Ro = 313m Q and the exponential rate coefficient, ap, = 
0.884 A~!. The values for the second double layer resistor are 
RR = 15.3mQ, RE = 199m Q and the rate coefficient, az = 
0.907 A~!. The fits are presented in Fig. 4(b). 

The EIS data also provides an indication of the values for the 
double layer capacitances and adsorption inductance. Literature 
indicates that these elements have constant values [16,13]. Fig. 5 
shows the values for the double layer capacitances (C)-mark for 
CpL and @-mark for Cp) and the adsorption inductance La 
(LJ-mark). 


Parameter values of the third-order equivalent circuit model of the PEM fuel cell based on the data from electrochemical impedance spectroscopy 


DC set-points 


1A 2A 3A 4A 5A 6A 
Roo(&) 0.283 0.270 0.290 0.257 0.266 0.265 
Cox (F) 0.221 0.259 0.263 0.243 0.229 0.226 
Rp_(mQ) 202 107 107 82.3 65.7 64.5 
Cr (F) 0.450 0.364 0.745 0.691 0.566 0.795 
Rr(m2) 92.3 54.9 18.1 21.7 24.2 18.2 
La (mH) 36.0 45.3 37.3 353 38.5 48.4 
Ra(mQ) 19.1 20.4 16.6 12.3 9.37 8.2 
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Fig. 4. Resistor values found with impedance spectroscopy and their fitted function results. In (a) a quadratic and in (b) a linear relationship between the fuel cell 


current and operating temperature is used. 


From Fig. 5, it is clear that the first double layer capacitor 
Cp. = 0.241 F and the adsorption inductor La = 40.2 mH can 
be interpreted as constants. The second double layer capacitor 
(CR), however, can hardly be assumed as constant. 

The second double layer capacitor is more difficult to fit in the 
higher DC current region, because the semi-circle in the complex 
impedance plane as the result of the parallel circuit Rp and Cr 
is dependent on the value of the second double layer resistor Rr. 
Because this resistor becomes small at higher output currents, 
the amplitude of the semi-circle also becomes small. 


3.2. Steady-state characteristic 


Steady-state measurement results are presented as C)-marks 
in Fig. 6. The characteristic shows almost no concentration 
polarization, which can be interpreted as the concentration coef- 
ficient that goes to zero, (B — 0 V) [10,9]. As a consequence, 
the expression for the fuel cell voltage (1) may be simplified to 


I , 
Epc = Eoy = A In (=) = R" Jpc. (9) 


The theoretical value of the reversible voltage of a single 
cell is 1.229 V under standard conditions (that is for tempera- 


C (F) , L (H) 


Fig. 5. Measured values and fitted as a constant (lines) of the two double layer 
capacitors (Cpy and Cpr), and the adsorption inductor (La). 


tures of 25°C and standard pressure of the reactants). However, 
because literature indicates that the standard operating temper- 
ature of a PEM fuel cell is around 60°C [8-10], and because 
the reversible cell voltage is dependent on the operating temper- 
ature [8,4,9,10], the reversible cell voltage is calculated to be 
Erev = 1.17 V. As was indicated in Section 2, the used fuel cell 
stack contains four membranes in series, resulting in a theoretical 
open circuit voltage of 


Eocv = 4E rey = 4.68. (10) 


With the values for the different resistors known, the steady- 
state resistance R'™' can be determined. Because in steady-state 
operation the adsorption inductance La is short circuiting the 
adsorption resistance Ra, the total internal resistance, which is 
dependent on the output current of the fuel cell stack (Jgc), is 
determined by 


R'™ pc) = Roo + Rpo-Urc) + RrUec). (11) 


Because there are two possibilities of determining the values 
for the double layer resistors, the results of the open circuit volt- 
age source, shown in Fig. 7, are also presented in two forms. 
Fig. 7(a) shows the open circuit voltage source (Epc) of the fuel 
cell stack for double layer resistors presented with (4) and (5), 


Fig. 6. Steady-state characteristic measured at the four membrane stack. The 
©-marks are the measured points, the full line is the fitted characteristic. 
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Epc (V) 


leg (A) 


Fig. 7. Open circuit source characteristics (O) and the fit (full line). In (a) a 
quadratic and in (b) a linear relationship between fuel cell current and operating 
temperature is used. 


while Fig. 7(b) presents Erc for the double layer resistors from 
(7) and (8). 

The difference between Fig. 7(a) and (b) directly shows the 
problem in fitting the double layer resistors with the current 
temperature relationship of (3). The open circuit voltage source 
of the stack has a minimum value at Ipc = 4A, while the Efc 
in Fig. 7(b) has a continuously decreasing trend. 

The large discrepancy between the first measurement point 
(Irc = 0A) and the simulated value, which can be seen in 
Fig. 7(b), is related to the internal current of the fuel cell [8]. This 
loss current prohibit the fuel cell from producing its open circuit 
voltage. A schematic representation of the fuel cell including 
the internal current source is shown in Fig. 8. 

In order to decide whether assumption (3) or (6) is best, 
the steady-state model (10) is first fitted using the least squares 
method to the open circuit source of Fig. 7(a). To include the 
minimum in the characteristic (at Ipc = 4A), the activation 
coefficient A has to be rewritten. The proposed open circuit 


Fig. 8. Representation of the fuel cell including the internal current source Jn. 


voltage source is 
I 
Epc = Eoey — (a + blĝc) In (=) (12) 
0 


where the activation coefficient is made dependent on the square 
of the fuel cell output current. The fit results in the constant 
term of activation a = 0.216 V, the activation resistive term b = 
541 uQ AT! and the exchange current Jọ = 972 pA. The error 
made between the open circuit function (the full line of Fig. 7(a)) 
and the measured data points is less than 2% except for Ipc = 
OA. 

Second, the kinetic mode steady-state model is fitted to the 
data from Fig. 7(b). For this fit, the activation coefficient A is 
assumed to be constant, resulting in 


Erc = Eocy — A In (=) : (13) 
Io 
where A = 0.164 V and Jọ = 117 p A. The error between this 
function (the full line of Fig. 7(b)) and the obtained data is less 
than 1%, except for Ipc = OA. 

Although both solutions for the double layer resistors are 
applicable to find a model for the open circuit voltage source 
Epc, solution (13) is more attractive. Not only is the error 
between the fitted open circuit voltage source in favor of the 
linear current dependency of the double layer resistors, but also 
the complexity of (13) is lower than that of (12). 


3.3. Step-response measurements 


With the step-response test, the load current of the fuel cell is 
suddenly changed from one value to another. This method can 
be useful in studying the fuel cell response to switched-mode 
converters, because the input of these converters would excite 
the fuel cell in a similar way. Step measurements fitted to a 
first-order model show a step polarity dependent time-constant, 
which is unsuited for a linear circuit model description [9,10]. 

By fitting the step-response to the proposed dynamic model 
of Fig. 9, including a constant value for Cp, the results indicate 
which value is to be used. The average value for all data points is 
Cr = 0.634 F (full line in Fig. 5), while the average value when 
omitting the last data point result in Cp = 0.562 F (dashed line 
in Fig. 5). 

For the verification of these assumed constant values for 
the capacitors and inductor, a simulation of the step-response 
is made and compared with the measured characteristic. In 
Fig. 10, the simulated and measured voltages of the fuel cell 
stack are shown in response to the current delivered. The simu- 


MN NVV NVV 
Ro | Rot | Rr 
Efc CoL Cr R 


- Ve + 


A La 


Fig. 9. Equivalent circuit model for a PEM fuel cell stack, featuring one or two 
double layer time-constants and an inductive adsorption time constant. 
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Fig. 10. Step-response: simulated voltage in Matlab (gray dashed line), mea- 
sured voltage (black line) and measured current (thin gray line). 


lated response is made while using Cr = 0.562 F for the second 
double layer capacitor. 

The error b etween the measured and simulated step-response 
is calculated using both modeling possibilities for the second 
double layer capacitor. The difference between the response for 
Cr = 0.562 F and Cr = 0.634 F is marginal, that is, less than 
3 mV peak and an average difference of less than 8 uV over the 
complete time interval. The largest voltage peaks of the error 
occur at the switching of the load. 

The step-response simulation is made by calculating the 
immediate changes in the resistances caused by the change in 
output current (2), (7) and (8) and combining those values with 
the steady-state voltage source (13), except for the value of the 
internal voltage source when Ipc = OA. For that specific situ- 
ation, an approximation of the internal voltage source is made, 
based on the value found in Fig. 7(b). 

The equivalent circuit model of Fig. 9 is determined with 
analytical expressions listed in Table 2. For the analytical equa- 
tions of the double layer resistors Rp, and Rr, and for the open 


circuit voltage source Epc, the linear relationship (6) between 
temperature changes and output current is used. 


4. Model verification 


The model, with parameters as given in Table 2, is verified 
with simulations in Matlab of the steady-state, small- and large- 
signal response, and the results are compared with the measured 
data. First, the electrochemical impedance data will be verified. 
Second, the steady-state characteristics are checked, and last, 
the step-response is inspected. 

The definition of the error (£) used in the verification is 


Vecsim — Vi 
—_ FC,sim FC,meas x 100%, (14) 


VFC,meas 


and is the relative error between the simulated voltage VfC,sim 
and the measured voltage VfFC,meas Of the fuel cell. 


4.1. Small-signal verification 


The simulation of the frequency response of the model is 
done by implementing the complex impedance in Matlab. After- 
wards, the calculated values of the impedance are compared to 
the measured ones. The relative error € between measurement 
and simulation is shown in Fig. 11. The top left graph shows 
the model error with a DC set point of Ipc = 1 A. The top right 
graph shows the error for Ipc = 2 A DC setting, the middle left 
for Ipc = 3 A, etc. 

From the figure, it can be seen that the largest deviation 
appears for | and 3A DC set point. This can be explained by 
referring to Fig. 2, where the same measurements show a slight 
shift of the semi-circles over the real axis. This shift results 
in a higher measured bulk resistance Roo than is used for the 
simulation, which results in a higher relative error. 


Table 2 
Analytical equations and values of the fuel cell model for a four membrane PEM fuel cell stack 
Model equations Components Values Units 
Erc = Eocv — A ln(Ipc/10) Eocv 4.68 V 
A 0.164 Vv 
Io 117 pA 
Roo Roo 0.272 Q 
CDL CoL 0.241 F 
RpL = RỌ + R? e7e2Ltrc RY 67.8 m2 
Ro 313 m2 
aDL 0.884 AT! 
CR CR 0.562 F 
Rg = RẸ + Re oriec R 15.3 m2 
RR 199 m2 
aR 0.907 AT! 
La La 40.2 mH 
Ra = RÌ — alec RY 23.5 m Q 
oA 2.62 mQAa! 
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Fig. 11. Relative error between the measured small-signal characteristic of the fuel cell and the simulated characteristic, for different DC set points. 


The error at Ipc = 4 A increases at higher frequencies. This 
effect is also related to a difference in bulk resistance compared 
to the used simulation value. An error in the bulk resistance 
will be enlarged in the higher frequency range, because the bulk 
resistor becomes a relatively larger part in the internal resistance 
of the model (R'"). This can be seen in the Ipc = 1 A, Ipc = 3 A 
and Ipc = 4A plots. 

From the results shown in Fig. 11, it can be concluded that 
the average error made in the small-signal behavior is lower than 
4%. The peak error is 6%. 


4.2. Steady-state verification 


The error (e) made in the steady-state simulation is pre- 
sented as the O-marks connected with the dotted line in Fig. 12. 
As predicted, the error between the measured output voltage 
at Ipc = 0A deviates most. The other verified points are very 
close to the obtained values, and the error differs at most 3% (at 
Ipc = 6A). 

To reduce the error at very low output currents of the stack, the 
steady-state equation (13) is expanded with an internal current 


0 1 2 3 4 5 6 


Fig. 12. Relative error between the measured and the simulated steady-state 
characteristic. 


term Jy, as proposed in reference [8]. The values for the open 
circuit voltage source (Epc) are fitted to 


Irc + =) 


T (5) 


Ec = Eoy — A In ( 
and for the new expression the parameters change to A = 
193 mV, Jo = 631 uA and the fuel crossover or internal cur- 
rent term will be J, = 470 mA. The fuel cell stack is made from 
four membranes with an active area of Ae = 60 cm? each, which 
leads to an internal current density of ja = 2mAcm~?. This 
value is marked in the literature as “normal” [8]. 

The results of the steady-state verification using the internal 
current lead to the error plot with the full line and the -marks 
of Fig. 12. The error is reduced to 3.6% maximum at Ipc = 6A, 
and the error at Ipc = 0A is almost 0%. 

The internal current component only influences the behavior 
of the steady state in the activation polarization region. In the 
higher current region, the influence of the internal impedance 
(R'") of the fuel cell becomes higher, while the influence of the 
activation polarization diminishes. 


4.3. Step-response verification 


The verification of the step-response is done in two ways. 
First, the error plot of the step-response presented in Fig. 10 
is simulated with Matlab with an open circuit voltage source 
without internal current term (J,). The top graph of Fig. 13 
shows the relative error as defined in (14) between measure- 
ment and simulation. The alternative, an open circuit voltage 
source (Epc) with Zn, is shown in the bottom graph of 
Fig. 13. 

The error of the step-response is highly dependent on the 
open circuit voltage source. As is seen in Fig. 13, the internal 
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Fig. 13. Relative error between the measured and simulated step characteristics. 
In the top figure, the FC voltage is simulated without internal current, while in 
the bottom figure it is simulated with internal current. 
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Fig. 14. Absolute error between the measured and simulated step characteristics. 


current term in the steady-state characteristic (15) influences the 
large-signal behavior, especially in the low output current region 
(activation polarization). It can also be seen in Fig. 13 that the 
errors of the large-signal behavior at Ipc = 6 A, that is the steps 
in the time interval 155 < t < 180, are higher when the internal 
current is included than without this term. This is also related to 
the error in the steady-state behavior, as can be seen in Fig. 12. 

The difference in absolute value between the simulated volt- 
age and the measured characteristic is very small over the 
complete current range, as is shown in Fig. 14. However, the 
error in the higher current range increases relatively more 
because of a lower fuel cell output voltage, as can be seen in 
Figs. 10 and 14. 

The large peak errors in the activation polarization region 
seen in Figs. 13 and 14, are caused due to instantaneous 
changing of the resistor values as a result of a different load 
current. In practice, the change in resistance will not be so 
abrupt. The peaks have a limited duration of only a couple of 
microseconds. 


5. Conclusions 


This paper provides a way of constructing an accurate fuel 
cell model based on electrochemical phenomena. The model is 
constructed as an equivalent circuit model, and is capable of 
simulating the static and dynamic behavior of PEM fuel cells, a 
so-called dynastatic model. 


The results obtained from electrochemical impedance spec- 
troscopy are used to find the small-signal behavior. From the EIS 
data the values of the resistor parameters can be found, together 
with a first indication for the capacitors and inductor values of 
the third-order equivalent circuit. The two double layer resistors 
(Rp and Rp) are non-linearly dependent on the output current 
(Irc). The bulk resistor (Roo) may be assumed constant and 
the adsorption resistor (Ra) is linearly dependent on the output 
current. 

Once the internal resistance is known, the open circuit volt- 
age source (Epc) of the equivalent circuit can be retrieved. In 
order to reduce the simulation error at no load, an estimation 
of the internal current (Zn) of the fuel cell is introduced. The 
additional internal current parameter is not only necessary to 
reduce the steady-state error, but also to reduce the error in the 
step-response. 

Finally, the indicative values of the capacitors and inductor 
from EIS can be used to find the optimal values with respect 
to the step-response measurements. When the capacitors and 
inductor values are fitted, the dynastatic model is complete. 

The fuel cell model is characterized by performing small- 
signal, steady-state and large-signal measurements. Because 
the modeling approach fully relies on these measurement 
results, electrochemical phenomena like temperature depen- 
dency, charge double layer and adsorption effects are implicitly 
covered with the current dependent element values. 
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